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Abstract 

We study the finite temperature crossovers in the vicinity of a zero tem- 
perature quantum phase transition. The universal crossover functions are 
observables of a continuum quantum field theory. Particular attention is fo- 
cussed on the high temperature limit of the continuum field theory, the so- 
called "quantum-critical" region. Basic features of crossovers are illustrated 
by a simple solvable model of dilute spinless fermions, and a partially solv- 
able model of dilute bosons. The low frequency relaxational behavior of the 
quantum-critical region is displayed in the solution of the transverse-field Ising 
model. The insights from these simple models lead to a fairly complete un- 
derstanding of the system of primary interest: the two-dimensional quantum 
rotor model, whose phase transition is expected to be in the same universality 
class as those in antiferromagnetic Heisenberg spin models. Recent work on 
the experimental implications of these results for the cuprate compounds is 
briefly reviewed. 
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I. INTRODUCTION 



Consider a quantum system on an infinite lattice described by the Hamiltonian Ti-lg), 
with g a dimensionless couphng constant. For any reasonable g, all observable properties 
of the ground state of 7i will vary smoothly as g is varied. However, there may be special 
points, like g = gc, where there is a non-analyticity in some property of the ground state: we 
identify gc as the position of a quantum phase transition. In finite lattices, non-analyticities 
can only occur at level crossings; the possibilities in infinite systems are richer as avoided 
level crossings can become sharp in the thermodynamic limit. In this paper, I will restrict 
my discussion to second order quantum transitions, or transitions in which the correlation 
length and correlation time diverge as g approaches gc- As I will review below, any such 
quantum transition can be used to define a continuum quantum field theory (CQFT): the 
CQFT has no intrinsic short-distance (or ultraviolet) cutoff. The main purpose of this paper 
is to review some recent wor on the properties of Ti.{g) at finite temperatures (T) in the 
vicinity of g = gc- This is equivalent to a study of the finite T crossovers of the associated 
CQFT. We shall focus especially on the dynamic properties of a ubiquitous finite T region, 
usually called "quantum critical"! (as we shall see below, there are reasons why this name is 
misleading and not quite appropriate; nevertheless, I will use it here) . The quantum-critical 
region appears as the high T limit of the CQFT; unlike the statics of classical lattice models, 
the high T limit of a CQFT is usually highly non-trivial. All of this discussion will take 
place in the context of some simple examples drawn from quantum spin systems. 

I set the stage by reviewing the Wilsonian approach to critical phenomena and field the- 
orieJ^, using the perspective of quantum critical phenomena. By the usual Trotter product 
decomposition, we can set up the partition function of Ti-lg) as a functional integral over 
degrees of freedom which fluctuate as a function of the spatial co-ordinate x and imaginary 
time T. Let us now examine the behavior of this functional integral under the rescaling 
transformationiH 

X e~^x T e^'^V (1) 

The dynamic exponent z determines the relative scaling dimensions of space and time co- 
ordinates. The critical point at g = gc is a fixed point of ([^), and g — gc is a relevant 
perturbation away from this point. We have therefore the flow equation 

which defines the critical exponent u- (For simplicity I do not discuss the case of fixed points 
with more than one relevant perturbation, as they can be treated in a similar manner). In 
the long-distance, long-time limit, this deviation from the critical point will be characterized 
by some renormalized energy scale, G. I emphasize that G is a dimensionful parameter, 
expressed in the laboratory units of energy, and directly measurable in an experiment; a 
typical example would be an energy gap. In the vicinity of the critical point, the renormalized 
energy scale G will be related to the bare coupling g by 

Gr^K\g-gcr (3) 
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where A is an ultraviolet cutoff, measured for convenience in the units of energy too. From 
the perspective of a field theorist, the CQFT associated with the quantum critical point is 
now defined by taking the limit A — >• cx) at fixed G; from (|]) we see that, because zv > 0, 
it is possible to take this limit by tuning the bare coupling g closer and closer to the critical 
point as A increases. (A condensed matter physicist would take the complementary, but 
equivalent, perspective of keeping A fixed but moving closer to criticality by lowering his 
probe frequency ~ G). The resulting CQFT then contains only the energy scale G. At 
finite temperatures, there is a second energy scale T (using units in which ks = 1); its 
thermodynamic properties will then be a universal function of the only dimensionless ratio 
available — G/T. It is the purpose of this paper to review recent work on the crossovers as 
a function of G/T in a number of systems, and to highlight the unusual properties of the 
heretofore unexplored high-temperature, "quantum critical" , limit of the CQFT, T 3> G. 

It is now easy to see why the high T limit of the CQFT can be non-trivial. A conventional 
high T expansions of the lattice model Ti proceeds with the series 

Tre-^/^ = Trl - ^TiH + + • • • (4) 

The successive terms in this series are well-defined and finite because of the ultraviolet 
cutoffs provided by the lattice. Further, the series is well-behaved provided T is larger than 
all other energy scales; in particular we need T ^ A. In contrast, the CQFT was defined 
by the limit A — > oo at fixed T, G, so the high T limit of the CQFT corresponds to the 
intermediate temperature range G <^ T <^ A of the lattice model. It is not possible to access 
this temperature range by an expansion as simple as (^), and more sophisticated techniques, 
to be discussed here, are necessary. (One could also, of course, determine a large number 
of terms in (^ and then use some Fade extrapolation methods to access the T ^ A region: 
this method has been used by Sokol et. alM and I will not discuss it here). 

In contrast to the static properties, the dynamic properties of Ti. are already non-trivial 
in the high T 3> A limit of the lattice model. Although one expects some sort of incoher- 
ent, dissipative dynamics, the damping co-efficients cannot be determined directly — all the 
approaches used so far are essentially variants of the methods discussed by Moriya0 and 
Forsterlli, and use a short-time expansion, coupled with an ansatz for the spectral function, 
to extrapolate to the long-time limit. In this paper, we will discuss the dynamics in the 
high T limit of the CQFT, G <C T -C A, or the "quantum-critical dynamics". The dy- 
namics continues to be dissipative and relaxational, but is not amenable to a description 
either by a classical Boltzmann equation for a dilute gas of quasiparticle excitations, or by a 
classical Langevin-like models of the types discussed in the classic review of Hohenberg and 
Halperinlli. However, as we shall show, the scaling structure of the CQFT does permit some 
progress to be made; indeed we will discuss below the complete solution of the quantum- 
critical dynamics of a simple spin model, including the exact determination of a damping 
co-efficient. 

We will begin our discussion with a simple solvable model of spinless fermions in Sec- 
tion |I|: this will allow introduction of the main concepts in a very simple setting. Section |T| 
will extend these results to a related but more complex model of dilute bosons in spatial 
dimensions d < 2. An explicit solution of the relaxational dynamics of the quantum-critical 
region will be obtained in the discussion in Section |rV| on the Ising model in a transverse 
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field. The expository analyses of these toy models will lead in Section to the main system 
of interest — the 0(3) quantum rotor model in o? = 2. We will also review applications of 
these results to numerical simulations and experiments. Finally Section |VI| will conclude by 
highlighting the main results and noting recent work on related subjects. 

II. DILUTE GAS OF SPINLESS FERMIONS 

Much of the physics I wish to discuss is displayed in a surprisingly simple model of a dilute 
gas of spinless fermions at finite temperature: its scaling forms have a structure identical to 
those of much more complicated models. The main shortcoming of the model is that the 
associated CQFT has no interactions, and there is therefore no relaxational behavior in the 
scaling functions. 

Consider the following Hamiltonian: 

1-Lf = -tY^ {c\c.j + c]ci - c\ci - c]cj) - (5) 

<ij> i 

where is a spinless fermion annihilation operator at the site i of a c?- dimensional hypercubic 
lattice, and < ij > are nearest neighbors. There are no interactions in Tip, so it is trivially 
solvable. Consider the ground state of Hp s.s a function of the dimensionless coupling 
constant g = fi/t. For g < gc = 0, the ground state has no particles. There is a non- 
analytic onset in the density of particles a.t g = gc, signaling a quantum phase transition. 
The bandwidth ~ t of the fermions plays the role of the upper cutoff in energy (A) for this 
transition, and the critical region defining the applicability of a CQFT is roughly T, /i -C t. 
In fact, it is not difficult to determine the exact effective action of the CQFT (in units with 
h = 1): 

Cp = J d^x dT^Ux,r) i — - — ^p{x,t) (6) 

where \E'ir(xj) = a~'^Ci, m ~ l/{ta?), and a is the lattice spacing. For this action we can 
identify by the usual methods the exponents v = 1/2 and 2; = 2 at the /i = critical point; 
taking the density n = {"^^p'^p) as the order parameter, we get the exponents f3 = d/2 and 
f] = d. Further, it is also easy to see that all interactions are irrelevant at this critical point 
in all dimensions d > (the least irrelevant interaction term, |\E'j7V\E'i?p, becomes relevant 
only for d < 0). Finally, for the renormalized energy scale measuring the deviation from the 
critical point, G, we may take G = n, the bare chemical potential in Tip. Note that there is 
no non-universal scale-factor in the relationship between G and fi — this is a consequence of 
the triviality of the critical exponents. More typical models with anomalous exponents will 
have non-universal scale-factors. 

We show in Fig |I|the phase diagram of T^i;' as a function of G = fi and T. Our interest is 
primarily in the regions within the hatched lines where the CQFT (^ applies. Within this 
region there are three physically distinct types of behavior (A, B and C): as Cp is trivially 
solvable, the universal properties of A, B, and C and the crossovers between them can be 
exactly determined. Let us describe the regions in turn: 
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(A) Activated fi ^ — T: The fermions are very dilute, with a density ~ e^^'^. Quantum 
effects are suppressed and the particles behave classically. 

(B) Fermi or Luttinger liquid /i ^ T: Now quantum degeneracy effects are paramount. At 
T = 0, the ground state is a Fermi liquid (in c? = 1, a Luttinger liquid); at finite T thermal 
effects lead to a small number of particle and hole excitations near the Fermi surface. 

(C) Quantum Critical <^ T: Unlike A and B, the temperature T is the most important 
energy scale in this region. We can set /x = here without much damage (all corrections 
involve positive powers of /i/T). The energy of a typical excitation in this region is of order 
T and as a result, quantum and thermal fluctuations are equally important. 

The relationships between the regions becomes clearer upon considering an explicit ex- 
ample of a crossover function. The density of particles n = (^^^j?) obeys the scaling 
form 

n = (2mT)^/2$. (^) (7) 
where the universal scaling function ^n{j^) is given by 

^niH) = J (^^^.^- (8) 

Notice that there are no arbitrary scale-factors in (0). In the activated region A (/I -C — 1), 
we have $ri ~ e'^/ (4??)'^/^, which is exactly the result we would have obtain from the classical 
Maxwell-Boltzmann statistics. In the Fermi liquid region B (/I ^ 1), $„ ^ Jt^^'^Vd/ {27iY, 
where is the volume of the unit sphere in d dimensions; this is the fully quantum result 
obtained by filling up the Fermi sphere. Most interesting is the quantum critical region C 
(|7i| <C 1) where 

M-p) = c[fj{l)\l-2'-'^')+OiJl) (9) 

The value of $„(0) depends upon the details of the Fermi distribution function, and not just 
its forms in the classical and quantum limits: this illustrates our assertion that quantum 
and classical effects are equally important in region C. 

It is also interesting to compare the behavior of the density in the universal region C 
with the true high temperature limit of the lattice model - region D. In other words, we are 
going vertically upwards in T at /i = in Fig |I[ It is easy to compute: 




where Ci is a non- universal constant (recall that in this model G = fi and A = t). Notice 
the difference between the universal high T limit of the CQFT (the first result) and the 
lattice high T limit. I hope that this example has illustrated the general principle, and in 
the remainder of the paper I will make no further reference to non-universal regions like D. 
It will be implicitly assumed that we are working with the universal continuum theory, and 
we will describe only regions like A, B and C. 
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Before closing our discussion on this deceptively simple model, we discuss the scaling form 
of observables as a function of momentum k and frequency u. Because of the absence of 
interactions, the single-particle Green's function is trivial; so we discuss the density-density 
correlator Xn{k, uj) which has a slightly more interesting structure. As the particles are free, 
Xn is of course given simply by the Lindhard function, which can be manipulated into the 
scaling form 

Uk.u>)^^^^,„(^-^,-,-j (11) 

where is a universal complex-valued function related to the Lindhard function. Notice 
that, like the scaling form (|^, there are again no arbitrary scale factors. Further, is 
well-defined at /i/T = 0, where it yields the dynamic susceptibility of the quantum-critical 
region C. We will see several other examples of scaling forms like (|^) and ([Tl|) in this paper, 
but the scaling functions will not be as simple as they are here. 



III. DILUTE BOSE GAS 



Now we consider the same density onset quantum transition considered in Sec |T|, but for 
the case of bosons. The discussion here is drawn from that of Sachdev, Senthil and Shankar§ 
to which the reader is referred for further details. 

Unlike the case for spinless fermions, it is no longer possible to ignore the interactions 
between the particles. We consider the properties of the following continuum model 

Cb = J d^x drvj/Ua;, r) f — - — - j ^b{x, t) 

+ \j d'^X j d'^x' j^J^ dT\mB{x,T)\^v{x-x')\^!B{x\T)\' (12) 



where b is a boson annihilation operator, and v{x) is a repulsive interaction of range ~ a. 
Like Cpi J^B has a quantum phase transition at /i = T = 0, and we will discuss its universal 
properties here. Because of the finite range of f , the universality only sets in at distances 
larger than a. A straightforward RG analysis^ of the vicinity of the quantum critical point 
shows that v flows into the = fixed point for d > 2. It turns out that v is actually 
dangerously irrelevant for d > 2: we do not wish to enter into a discussion of such effects 
here, and so most of our remaining discussion will be restricted to d < 2. 

For d < 2, v{x) flows into a universal fixed point interaction v{x) = u*6{x) for some d 
dependent constant u*. The scaling structure of this fixed point turns out to be very closely 
related to the non-interacting spinless fermion model of Sec |I| for the same value of d. All 
exponents and scaling forms of the boson and fermion models are identical, but the scaling 
functions themselves are different. The crossover phase diagram of the d < 2 dilute Bose gas 
is essentially identical to the fermion phase diagram in Fig|l|, but the physical interpretation 
of the phases is somewhat different: 

(A) Activated fi ^ —T: This is essentially identical to the fermion case as the particles are 
dilute and their quantum statistics plays a negligible role. 
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(B) Incipient Superfluid fi ^ T: The ground state is now a superfluid (in (i = 1 a Luttinger 
liquid), but classical thermal fluctuations destroy the long range order at any non-zero 
temperature. Nevertheless, the phase coherence length is large and system behaves like a 
superfluid at short scales. 

(C) Quantum Critical ^ T: This is similar to the fermion case in that T is the most 
important energy scale. However, there are now strong interactions among the particles, 
leading to an incoherent excitation spectrum. The system does not display characteristics 
of a superfluid ground state at any length scale, but instead crosses over directly from free 
particle behavior at short time scales, to dissipative, relaxational dynamics at long time 
scales. The single particle Green's function G{x,t) = (\l'B(a;, r)\E'^(0, 0)) has a non-trivial 
scaling function $g (this form actually holds in all three regions A, B, and C) 



We have Fourier transformed and analytically continued to the retarded Green's function 
at real frequencies. This scaling form also holds for the spinless fermions of Section O, but 
the scaling function then is simply the free fermion form $0(^5 oJ, /l) = l/(a7 — — /I + irf) 
where 77 is a positive inflnitesimal. Computing $0, and other scaling functions, for bosonic 
quantum critical dynamics is not as easy. One approach! is to expand in powers of the 
flxed-point interaction u* which becomes small as d approaches 2 from below: this becomes 
an expansion 'm. e = 2 — d. 

In d = 1, it is possible to make more explicit progress. It has been argued! that now 
u* = 00. The bosons thence become impenetrable, and their quantum mechanics becomes 
identical to those of free fermions. Hence, in d = 1, the quantum critical dynamics of dilute 
gases of spinless fermions and boson are described by the same CQFT, Cp of Eqn (|]). For the 
bosonic system we have to supplement Cp with the following non-local relationship between 
the boson and fermion operators (essentially a continuum Jordan- Wigner transformation) 

"^six) = exp (^Tc dx'^Ux')^F{x')j'^F{x) (14) 

We are not home yet, as evaluating correlators of (|1^) under (|^) is not easy. Korepin and 
SlavnoAiiB have succeeded in showing how this problem may be reduced to determining 
the solution and Fredholm determinant of a linear Fredholm integral equation. At this 
stage, numerical analysis is required, and some scaling functions have been determined to 
essentially arbitrary accuracyS. 

Finally, an additional comment about the d >2 case. There is now true superfluidity and 
a flnite temperature phase transition to a normal state, all within region B. The remainder 
of the phase diagram remains the same as in the d <2 case. 



IV. ISING MODEL IN A TRANSVERSE FIELD 

Unlike the dilute Fermi and Bose gases, the Ising model possesses anomalous exponents. 
Yet it is simple enough in c? = 1 to allow exact computation of a quantum-critical dynamic 
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correlation function. We will also review, in this section, the very useful and general mapping 
between the quantum model and an equivalent classical statistical mechanics model; we will 
then discuss the crossovers in the phase diagram like Fig |I| in the context of the classical 
model. Some of of the following discussion is a review of well-known properties of the Ising 
mo our main purpose here is to use the explicit solution in ci = 1 to present a 

physical interpretation which generalizes to other quantum phase transitions. 

We will explicitly discuss the following Hamiltonian, describing the Ising model in a 
transverse field in c? = 1 (we will remark briefly on the generalization to higher d): 

Hi = -JY^ {ga^^i + a^^icr^^i+i) (15) 

i 

where J > in an overall energy scale, > is a dimensionless coupling constant, and 
c"a;,i, CTz,! are Pauli matrices on a chain of sites, i. Consider the ground states of Tii for small 
and large g in turn. For small g, the second term in (|T5|) dominates and the spins all align 
themselves either in the +z or —z directions: there is a spontaneous magnetization and spin- 
reversal symmetry is broken. On the other hand, for large g, the first term in ( [151 ) prefers 
a state which is in a superposition of eigenstates, with different sites uncorrelated: the 
wavefunction looks like nj(|+)i+ \ These two limits are separated by a phase transition 
g = gc (in fact, gc = 1, exactly, because of a self-duality property of (p!5|)). In this section, 
we shall discuss the finite T dynamic properties in the vicinity of g = gc- 



To begin, we recall the well-known facttlHis that the T = correlators of Tij are similar 
to those in the classical, two-dimensional, Ising model given by the partition function Tr 
with 

F = -KY^ a^^ia^j (16) 

<ij> 

where the sites i,j now lie on a square lattice (say). This classical model has a phase 
transition ai K = Kc- At the level of the critica continuum theories, the mapping between 
the classical and quantum models become exact^ . There is a simple relationship between the 
two-dimensional field theory, with classical degrees of freedom, describing F in the vicinity 
of = i^c and the one-dimensional quantum field theory describing Tij near g = gc- one 
simply identifies one of the spatial directions of the classical theory as an imaginary time, and 
then analytically continues the correlators to real time, to obtain observables of the quantum 
theory. This mapping leads immediately to some useful information. As the classical theory 
is spatially isotropic, the quantum theory has dynamic exponent z = 1. Further, in the 
classical model the (Jz correlator behaves like ~ in momentum spaceEl {p is the two- 
dimensional momentum of the classical model); analytically continuing this to real time, we 
obtain for the dynamic susceptibility x(^) ^) (this is the az Oz correlator of the quantum 
7i/ and k is now a one-dimensional spatial momentum) as a function oi d= \ momentum k 
and frequency oj: 

X(fc>^)= (^2^^^2)7/8 T = ^,g = gc. (17) 

Here Z and c (an excitation velocity) are non-universal constants. We plot Imx(A;, ti;)/ti; in 
Fig 0. Notice that there are no delta functions in the spectral density, indicating the absence 
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of any well-defined quasiparticles. Instead, we have a critical continuum of excitations. We 
can also compute the dynamic local susceptibility, 



J Zn 



(18) 



a quantity often measured in neutron scattering experiments: 



2r(7/8)r(5/8)c \ \uj 



) 



3/4 



T = 0,g = gc- 



(19) 



This quantity is the density of states of local spin-flip excitations, and has a divergence as 

Our discussion so far has been at T = 0, and let us turn now to non-zero T. A finite T 
translates into a finite size Lr = 1/T for the classical model along imaginary time direction. 
Periodic boundary conditions are imposed along the finite direction, and so the classical 
model has the geometry of a cylinder with circumference = 1/T. We now discuss the 
crossovers at finite T in the context of the classical model. At L,- = oo (T = 0) we can 
characterize the deviations from criticality by the correlation length C, ^ \g — gd^^ (for the 
quantum model we may choose our renormalized energy scale G = 1/^ ^ \g — gd). At 
distances shorter than ^, the spins display critical correlations characteristic of the point 
K = Kc] it is only at distances larger than ^ that they become sensitive to the value of 
K — Kc and display characteristics of the ordered (for K > Kc) or the paramagnetic (for 
K < Kc) phase. Now consider the effect of a finite Lr] there are two distinct possibilities: 
(i) < Lt- (for the quantum model, T < G): Moving from the shortest to largest length 
scales, the crossover from critical to non-critical (either ordered or paramagnetic) behavior 
still occurs at a scale ~ C,', the length scale Lr has little effect at this point. The effects of 
Lr only become apparent at larger scales, at which point it is permissible to use an effective 
model which characterizes the non-critical ground state. 

(a) Lr < ^ (for the quantum model, T > G): Now the short-distance critical fluctuations 
see a flnite size Lr before they have had a chance to become sensitive to — Kc- These 
critical fluctuations are quenched by flnite size effects in a universal way. The resulting 
non-critical theory then responds only weakly at the scale ^. Note that the system does not 
display characteristics of the ordered or the paramagnetic state, of the = oo system, at 
any length scale. 

The above arguments are summarized in Figs ^ and ^; notice that Fig ^ is quite similar 
to Fig |I|. In all three regions of Fig ^, at the largest frequencies, u, Hj displays the critical 
correlations of the g = gc point as depicted in Fig ^ (we have returned now to the language 
of the quantum model Hi). 

In regions A and B there is a crossover from these critical fluctuations, at an energy 
scale G, to the behavior of the ordered {g < gc) or paramagnetic {g > gc) ground state of 
Hi (see Fig Both ground states have a gap, and thermal fluctuations will lead to dilute 
gas of quasiparticle excitations. We expect that an effective classical model (like Glaubeiill 
dynamics, which is similar in spirit to the Langevin models of Hohenberg and HalperinEl) 
will provide a suitable description of these thermal fluctuations. In c? = 1, on the ordered 
side {g < gc), these quasiparticle excitations are the 'kink' and 'anti-kink' solitons; even 
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an infinitesimal concentration of tfiese is sufficient to destroy long-range order at any finite 
temperature. However, for d > 1, long-range order is not immediately destroyed: as a result 
there is finite temperature phase transition within region A where the magnetic moment 
disappears; this transition will be in the universality class of the (i-dimensional classical 
Ising mo del0. 

In the quantum-critical region C, the critical fluctuations are quenched by thermal effects 
at the energy scale T (see Fig Q an early analysis of finite T crossovers in the transverse- 
field Ising model by SuzukitZl failed to identify region C). This quenching is completely 
universal and will be described explicitly below. The system has had no chance to display any 
characteristic of either non-critical ground state at any frequency scale. At low frequencies, 
the system realizes a new quantum relaxational regime. It is this regime which is really 
characteristic of the region C, and not the high frequency critical behavior which is present 
in all three regions. It is in this sense that the name "quantum-critical" of region C is a 
misnomer. 

As in the case of the dilute Fermi and Bose gases considered earlier, the dynamic sus- 
ceptibility x{k, will satisfy a universal scaling form over the regions of Fig ^ 



X{k,u) = -—^^ {—,-,-] (20) 




In the following we will determine the leading term in in the quantum-critical region C 
i.e. we will present an exact expression for $^(/c, tJ, 0). 

In the classical model F at K = Kc and Lr = oo, we know from {\iT\) that the Green's 
function G{x,t) = {az{0,0)az{x,T)) ~ (x^ + t'^c^)~^^^; we are using the label r for spatial 
direction corresponding to imaginary time. We can now use a remarkable result of Cardyil, 
which relies on the conformal invariance of this critical theory, to obtain an exact result for 
G in a system with a finite L^: 

r(i/8)z 1 y^v 1 y^^ 

~ 2i3/87r3/4r(7/8)c \u) l^cosh(27rx/L,c) - cos(27rr/L,) J ^^^^ 

This result has been asserted earlier by an inspection of the partial differential equation 
satisfied by G^. Although results like (pH) have been known to conformal field theorists for 
some time, they usually interchange the roles of x and t i.e. they consider systems of finite 
spatial length L^, and infinite temporal length, so that the system is in its ground state. 
For us, the spatial extent is infinite, and there is a finite length along the r direction, 
with the correlator (^) periodic in r with period (such a perspective has also been 
discussed by Shankaiu and by Korepin et. a/.El). Note that it doesn't really make sense to 
talk about the long imaginary time limit, r ^ L,-. However, after analytic continuation to 
real time, the long time limit of the quantum problem, t ^ 1/T or -C T, is eminently 
sensible, and is precisely the new quantum relaxational regime that we wish to access. The 
analytic continuation is a little more convenient in Fourier space: we Fourier transform (^Tj) 
to obtain G{k,uJn) at the Matsubara frequencies Un and then analytically continue to real 
frequencies (there are some interesting subtleties in the Fourier transform to G{k,Un) and 
its analytic structure in the complex u plane, which are discussed elsewhereS). This gives 
us the universal function $^ in the quantum-critical region: 
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r — + i r — + i 

1 \1Q Att J \16 Att 

(47r)V4 /i5 zj+k\ ^ /15 uJ -F 

r — h z r — h I 

V16 Att V16 An 



(22) 



We show a plot of Im$^/cJ in Fig |^. This result is the finite T version of Fig Notice that 
the sharp features of Fig ^ have been smoothed out on the scale T, and there is non-zero 
absorption at all frequencies. We can also observe the crossover as a function of frequency 
claimed earlier in Fig Notice that for aJ, :§> 1 there is a well-defined peak in Im$^/cJ 
(Fig ^) rather like the T = critical behavior of Fig 0. However, for cJ, A; ^ 1 we cross- 
over to the quantum relaxational regime and the spectral density Im$^/cJ is similar to a 
Lorentzian around a; = 0. This relaxational behavior can be characterized by a relaxation 
rate F/j defined a 



d\nxiO,uj) 



duo 



(this is motivated by the phenomenological relaxational form x(0,Ci;) 
0(cj2))). From (|2D|) and (^ we determine: 

ksT 



R 



n 

2 tan — 
16 



(23) 

Xo/(l - iuj/Vn + 
(24) 



where we have returned to physical units. The ease with which this result was obtained belies 
(I claim) its remarkable nature. Notice that we are working in a closed Hamiltonian system, 
evolving unitarily in time with the operator e~*^^*, from an initial density matrix given by 
the Gibbs ensemble at a temperature T. Yet, we have obtained relaxational behavior at 
low frequencies, and determined an exact value for a dissipation constant. Such behavior is 
more typically obtained in phenomenological models which couple the system to an external 
heat bath and postulate an equation of motion of the Langevin type. Notice also that Tii 
in ([T5|) is known to be integrable with an infinite number of conservation lawsS. However, 
the conservation laws are associated with a mapping to a free fermion model and are highly 
non-local in our degrees of freedom; they play essentially no role in our considerations, 
and do not preclude relaxational behavior in the variables. 

For completeness we also present results on a related observable which shows the crossover 
from critical to quantum relaxation behavior. We consider the local susceptiblity Xl? and 
obtain the finite T form of (p!9| ) by integrating 



21) over momenta: 



cT3/4 



UJ 

r 



T 



$l(cJ,0) 
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■ sinh ( — 
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(25) 
(26) 



27/%5/4r(5/8)r(7/8) 

A plot of the scaling function $l is shown in Fig ^ The function has the asymptotic limits: 

r^(i/8) 



$l(^,0) 



2iV%5/4r(5/8)r(7/8 



UJ 



sgn(a;) 



2r(7/8)r(5/8 




\uj\ < 1 



IcjI > 1 



(27) 
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The small a; behavior is relaxational as $l is linear in frequency, and the critical behavior 
at large frequencies agrees with ([T9|) . 

V. QUANTUM ROTORS IN TWO DIMENSIONS 

The study of this model is of direct experimental interest, as it is believedi to be a 
reasonable model of the spin fluctuations in antiferromagnetic compounds like La2CuOi 
and its lightly doped variants. The insight gained from the simple models studied in the 
previous sections will now be of great use, and we will rapidly be able to present a scaling 
analysis of its quantum phase transition. 

The Hamiltonian of the quantum rotor model is 

nn = JY.^-Ll-JY.na, (28) 

i <ij> 

where J > is an overall energy scale, (7 > is a dimensionless coupling constant, and i,j 
are the sites of a two dimensional lattice (< ij > denotes nearest neighbors). Notice the 
similarity between the forms of Hr and TCi in (|15|): it will turn out that the corresponding 
terms play a similar role. On each site i of the lattice we have the 3-component vector 
operators L, n (dropping the site index), which obey the commutation relations: 

[na,nb]=Q , [La,nb]=ieabcnc , [La,U] = ieabcLc- (29) 

The vector n is of unit length = 1, and its orientation identifies direction of the local 
magnetic order; the quantum rotor model is usually considered as an effective model for an 
underlying system of Heisenberg spins — in this case the magnetic order can be any ordering 
which is specified by a single vector and has no spatially averaged magnetic moment. The 
simplest example of this is the two sublattice Neel ordering, and we will therefore refer to 
(n) as the Neel order parameter. The L operator measures the angular momentum, and as 
all phases have no net magnetic moment, we will always have (L) = 0. 

For further insight into the meaning of Tin, consider the eigenstates of a single site 
Hamiltonian JgL'^/2. This describes a particle moving on a unit sphere with angular co- 
ordinate n and kinetic energy JgL'^/2. Its eigenenergies are Jgi{i + l)/2 with degeneracy 
2£ + 1 where i = 0, 1, 2, 3 . . .. The ground state is a non-degenerate singlet {i = 0) and 
has maximum uncertainty in the orientation of n. For large g, the ground state of Tin can 
be approximated by the tensor product of £ = states on each site. This state is clearly a 
quantum paramagnet and has a gap. A, to all excitations. Notice the similarity between this 
state and the large g quantum paramagnet of the Ising model Hi] in both cases the order 
parameters a^, n are in a state of maximum uncertainty. The small g limit of is also 
similar to the small g limit of Tij: now the exchange interactions between the sites prefer a 
state in which n has the same definite orientation on each site. Therefore, we expect long- 
range Neel order in the small g ground state. These two limiting states will be separated 
by a quantum phase transition at g = gc, which is, of course, the main subject of interest in 
this section. 

Before discussing the critical properties, we pause to remark on the relationship between 
the rotor model and Heisenberg antiferromagnets. Consider a pair of antiferromagnetically 
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coupled spin-S' Heisenberg spins: the eigenstates of this pair will have energies oc i{i + 1) — 
const for i = 0,1,... 25". Notice the similarity between these states and those of a single 
quantum rotor; the only difference is that there is no upper limit on the maximum value 
of i in the rotor case. However, it is reasonable to expect that these extra high energy 
states will not modify the low energy properties of lattice models. Therefore, there is little 
reason to doubt that the critical properties of Heisenberg antiferromagnets with a natural 
pairing of spins (or more generally, a natural clustering into an even number of spins) will 
be same as those of the rotor model. Antiferromagnets with no such pairing contain net 
Berry phase terms in their imaginary time path integral, beyond those present for the rotor 
model. However these Berry phases cancel between the sites, except for "hedgehog" -like 
spacetime singularitiesil, and it has been argued0i that these remnant Berry phases have 
no effect of the leading critical singularities. The reader is referred to the original papers for 
further discussion on these subtle issuesili: we will restrict our discussion here to the much 
simpler rotor model 7Y/j. 

The critical properties and phase diagram of the d = 2 rotor model Tin turn out to be 
remarkably similar to those of the transverse field Ising model Tij of ( [T5| ) in d = 1. Like the 
d = 1 Ising model, the d = 2 rotor model has no phase transition at any finite T: so the phase 
diagram of Fig |^ applies to 'Hr., with no phase boundary in region A (the phase diagram for 
Tiji m d = 2 was obtained first by Chakravarty et. al%). In both systems, the T = 0, g = gc 
critical point has z = 1. In the case of Tij this critical point was described by a CQFT which 
upon analytic continuation to imaginary time was the field theory of the two-dimensional 
classical Ising model. The analogous mapping for Tin yields the CQFT associated with 
the field theory for the three-dimensional, classical, Heisenberg ferromagnet. The universal 
scaling functions describing the crossovers in Fig |^ and ^ have an identical form in both 
theories, although, because the critical field theories are different, the critical exponents, 
universal amplitude ratios, and scaling functions will have different numerical values. The 
explicit results presented in Sec |V| for quantum-critical region C of the Ising model, all 
apply, unchanged in form, to the region C of the d = 2 rotor model: the qualitative features 
of the spectral functions in Figs ^, ^, and |^ remain the same, and the relaxation rate F^ 
(defined in (|2^)) satisfies ( p^ but with a different universal numerical prefactor. However, 
unlike the d = 1 Ising model, we cannot now get exact numerical results for the scaling 
functions of Ti.R: this is because the three-dimensional classical Hiesenberg ferromagnet is 
not exactly solvable (unlike the two-dimensional classical Ising model). Instead we have 
to be satisfied by approximate methods; reasonably accurate numerical estimates can be 
obtained in the 1/A^ expansion which has been discussed at length by Chubukov, Sachdev 
and Yei. 

There is a small, but significant, difference between the Ising model ind = 1 and the rotor 
model in d = 2 which cannot go unmentioned. This difference applies mainly to "ordered" 
regime in region A (See Figs ^ and H). The T = ground state of Hi for g < gc has a gap, 
associate with the finite energy cost of creating a kink or anti-kink soliton. In contrast, the 
g < gc ground state of T-Cr has gapless spin- wave excitations because of the broken continuous 
0(3) symmetry of the ordered state. So we can no longer use the gap. A, as the energy 
scale, G = A {gc — gY^ for measuring deviations from g = gc for g < gc- A convenient 
substitute turns out to be the spin stiffness G = ps which has the physical dimensions of 
energy in d = 2, and which also vanishes as Ps ~ {gc — g)^'^- At finite T in region A, it 
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is knownizli that the spin correlation length ~ exp{2'n'ps/T) as T — > 0. It is interesting to 
note that the behavior of Tij in c? = 1 is very similar: in this case the correlation length is 
determined by the mean spacing between kinks, and therefore behaves as ~ exp(A/T) for 
low T in region A. (Related to the exponential divergence of the correlation length, there is 
a further sub-division of the "ordered" regime of region A (Figs ^ and ^) at the energy scale 
c/ (correlation length); this complication occurs for both the d = 2 rotor and the d = 1 Ising 
models, and has been discussed elsewherei'i.) 

Associated with the continuous 0(3) symmetry of Ti/j, there is an important observable 
whose properties cannot be deduced by an analogy with the the Ising model. This is the 
uniform susceptibility xh, the response to a field, H which couples to the global conserved 
charge associated with the continuous symmetry: 



(30) 



The scaling dimension of xh can be determined exactly using symmetry arguments and the 
assumption of hyperscalingpa: this yields the scaling form 



T 



Xh 



H 



(31) 



Here c, is the same velocity that appears in a T = correlator like {\17\) , and is a fully 
universal function. This function was computed exactlyi'i in a 0{N = oo) rotor model, 
along with 1/A^ corrections in some limits; it has the limiting behavior 




0.6189 

— + 



+ . . . r ^ 



(32) 



r —>■ oo 



The two terms in the second result (r oo) are expected to be exact to all orders in 1/A^; 
the same two terms were also obtained by Hasenfratz and Niedermayer^l. Notice from (|3lD 
and ( p2D that xh has a linear dependence on T both for T <^ ps and T ^ ps] the slopes 
however differ by a factor of about 3 (for N = 3) and this will be important for experimental 
comparison 



A. Comparison with simulations and experiments on Heisenberg antiferromagnets 

The most straightforward comparison is with the double-layer, spin-1/2 Heisenberg an- 
tiferromagnet@. This model consists of spin-1/2 Heisenberg spins on two adjacent square 
lattices, with an intralayer antiferromagnetic exchange J and an interlayer antiferromagnetic 
exchange K. The ratio K/ J acts much like the dimensionless coupling g, with the large K/ J 
a gapped quantum paramagnet of singlet pairs of spins in opposite layers, and the small K/ J 
magnetically ordered. Extensive numerical simulations have been carried out on this model 
by Sandvik and collaborators^, and the critical point K = identified rather precisely. 
It is then possible to study the quantum-critical region C quite carefully as it extends over 
the maximum T range. A number of universal amplitude ratios, including those associated 
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with xhi and all results are now in good agreement with the expansion on the 0{N) 
quantum rotor model. Results from high temperature series expansions on the double-layer 
model also support this conclusion^. 

Secondly, comparisons have been made with the single-layer, square lattice spin-1/2 
Heisenberg antiferromagnet, both via simulations and by experimental measurements on 
La2Cu04,. This model has long-range order at T = 0, so must map onto the non-linear 
sigma model with g < g^. The low T region A was studied in the paper of Chakravarty, 
Halperin and Nelsoni, with good experimental agreement. Here we focus on the issue of 
whether this lattice model exhibits the CQFT high T behavior of region C, or it goes directly 
from region A to a non-universal, lattice dominated high T region like D of Fig ^. It was 
first argued by Chubukov and Sachdevi that this model does indeed possess a significant 
intermediate temperature regime of region C: this was based on comparisons with the T 
dependence of xhi with the factor of 3 alluded to above playing an important role. They 
also noted that it would be difficult to identify this region in the correlation length, an 
observation that was subsequently re-iterated by Greven et. a/.ii. A rather convincing 
demonstration of the presence of region C was given recently by Elstner et. a/.0'0 who 
examined a large number of observables in a high T expansion, and found good consistency 
with the universal rotor model results. 

Also significant in this context have been nuclear magnetic resonance experiments of Imai 
et. a/.iion La2Cu04. They have measured the T dependence of the 1/Ti and I/T2 relaxation 
rates at intermediate temperatures. Their observations are in reasonable agreement with 
the predictionliSii that can be obtained from the universal scaling result for x{k,^^) in 
the quantum-critical region C. 



VI. CONCLUSIONS 



This paper has presented a discussion of the vicinity of a second-order quantum phase 
transition in the context of a number of simple models. The overall picture that emerges 
is summarized in Fig |^, which shows a generic phase diagram in the plane of a coupling 
constant g, and the temperature T of a d-dimensional system; this phase diagram is a 
generalization of a diagram obtained first by Clralcravarty et. a3 for tire d = 2 quantum 
rotor model. The quantum phase transition occurs at the point g = gc, T = 0. Associated 
with this critical point, we can define a continuum quantum field theory (CQFT) over 
spacetime. In general, space (x) and time (r) do not play a similar role in the CQFT and 
have different scaling dimensions: x ^ x/s, r —>■ t/s^ under a spatial rescaling by s, with 
z the dynamic exponent. Further, even in imaginary time, the action for the CQFT can 
be complex due to the presence of Berry phases, and therefore corresponds to a statistical 
mechanics model with complex weights. Correlators of the CQFT are the universal functions 
describing crossovers in the vicinity oi g = g^ T = 0. The CQFT has no ultraviolet cutoff, 
but is characterized solely by two energy scales: the temperature T and an energy scale 
G ~ \g — QcY'^ characterizing the deviation of the ground state from the critical point (here 
V is the correlation length exponent). The value of the ratio G/T determines two distinct 
regions of the CQFT shown in Fig |^. In both regions there is a high frequency regime 
{uj > max(T, G)) which is dominated by excitations of the critical CQFT of the g = gc, 
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T = point. The regions are distinguished only by their low frequency behavior, which we 
discuss in turn: 

(z) The low T region (T ^ G, shown shaded in Fig |^) is a region of "conventional" physics 
for frequency scales uo < G. This frequency regime can be understood by beginning with 
the non-critical ground state and examining the particle-like excitations above it. A simple 
classical model (a Boltzmann equation for a gas of quasiparticle excitations, or a Langevin 
model of the types discussed by Hohenberg and Halperinll^l) is usually adequate for describing 
the long-distance, long-time dynamics of these excitations. The shaded region can also 
contain thermally driven, phase transitions; these transitions will be described by a classical 
field theory. 

(m) The high T region of the CQFT has a novel quantum relaxational regime. This regime 
is not described by an effective classical model, and displays intrinsic quantum-mechanical 
effects at the longest time and distance scales. A scaling analysis for this regime was reviewed 
in this paper. In some cases, as in the model of Section it is possible to obtain an exact 
value for the relaxation constant. 

This paper has reviewed only a small portion of what is a rapidly developing subject. 
We list below a number of recent (and not so recent) developments in related areas: 

• Quantum transitions between Fermi liquids and states with various types of spin or 
charge density wave orderings were discussed in important early work by HertJ. Hefo- 
cussed on the immediate vicinity of the finite termperature transition, like that within 
region A in Fig ^. In particular. Hertz missed the existence of the "quantum-critical" 
regime (as was pointed out recently by MillisH — this oversight is similar to Suzuki's0 
for the Ising model), which has been the main focus of this paper. More detailed 
studies of quantum transitions involving Fermi liquids have appeared recentlyS'&^i. 
A related, but different, perspective is provided by studies of critical phenomena in 
rotor models with doped electronsiS. 

• Related ideas on scaling in the quantum critical region have been presented by Tsvelik 
and collaborators!^. 

• A great deal of work has been done recently on "quantum impurity" models^El like 
the multi-channel Kondo effect. These models also display quantum phase transitions, 
with crossovers bearing some similarity to those discussed here. However the transi- 
tions do not modify the bulk properties, and are related instead to boundary critical 
phenomena. 

• As we indicated briefly in the discussion on the Bose gas in Section |ITT|, dangerously 
irrelevant operators sometimes need to be considered, as they do in the dilute Bose 
gas for d > 2. In fact such effects arise somewhat more frequently than they do in 
classical critical phenomena, as the upper critical dimension of the quantum transition 
is often quite low. An early analysis of the dilute Bose gas in c? = 3 by Weichmann 
et. a/.0 was dominated by such effects, although they did not present their results in 
the general context of quantum phase transitions. Such a perspective can be found in 
more recent workilS. 
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An important subject, on which much is not understood, is the effect of quenched 
randomness on quantum phase transitions. An early analysis for the random quantum 
rotor model was given by Boyanovsky and Cardyi. A recent exact solution by FisheiS 
of the random transverse-field Ising model in c? = 1 represents significant progress. The 
nature of spin glass ordering at T = and its destruction by quantum fluctuations has 
also been studied recentlySil. 

A number of experiment J^i^ have reported scaling of the type in Eqn in the local 
dynamic susceptibility. However the universality classes controlling these systems are 
not understood and quenched randomness appears to play a signiflcant role. Never- 
theless, it is interesting to note the qualitative similarity between the experimental 
measurements in Fig. 4 of Aronson et. a/.ii and our result for the local susceptibil- 
ity in Fig ^. The latter measures spin correlations on a single Ising spin induced via 
its coupling to its environment of other Ising spins. However, there is a fundamental 
equivalence between the spin being measured and its environment; this is an important 
difference between our bulk approach and alternative descriptions of the experiments 
using "quantum impurity" modelsSSii which clearly distinguish between the impu- 
rity and environment degrees of freedom. 

A recent studyil has examined the CQFT of a quantum ferromagnet. This system does 
not display a quantum phase transition of the type discussed here. Nevertheless, its 
phase diagram has regions similar to those in Fig ^ and aspects of its scaling properties 
are related to those of the dilute Bose gas of Section PTT| . 
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FIG. 1. Phase diagram of the dilute Fermi gas Hp (Eqn (^)) as a function of the chemical 
potential [i and the temperature T . The regions A, B, C, are separated by crossovers; all observables 
in these regions are described universal crossover functions of the CQFT Lp (Eqn (^)). Region A 
has an exponentially activated fermion density, region B is a Fermi liquid (in c? = 1, a Luttinger 
liquid), and region C is quantum critical. The hatched region marks the boundary of applicability 
of the CQFT and occurs at T ~ A = t. The same phase diagram also applies to the dilute Bose 



gas of Sec. ID, but the interpretation of the regions is different. 
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FIG. 3. Phase diagram of the transverse field Ising model in d > 1 dimensions; dashed lines 
denote crossover boundaries, while the solid line is a phase transition. In region A the system 
displays characteristics of the ordered ground state; thermal fluctuations about this state can 
destroy long-range order at any non-zero T in d = 1, but for d > 1 there is a phase transition 
(in the universality class of the d-dimensional classical Ising model). In region B the system is a 
gapped quantum paramagnet, and region C is quantum-critical. The crossover boundaries obey 
^ ~ Id ~ 9c\^'^ with z = 1 and i' the correlation length exponent of the d+ 1 dimensional classical 
Ising model (z^ = 1 in d = 1). The same phase diagram also applies to the quantum rotor model 
of Section ^ in d > 2, with the difference that the phase transition in region A is present only for 
d > 2 and is then in the universality class of the d+1 dimensional classical Heisenberg ferromagnet. 
Chakravarty et. a3 obtained this phase diagram for the d = 2 quantum rotor model. 
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FIG. 4. Crossovers as a function of probe frequency in the regions of Fig. ^ for the trans- 
verse-field Ising model in d dimensions. The d+1 critical regime is described by the CQFT at 
g = Qc and T = 0. The "ordered" regime is where a classical description in terms of a gas of 
well-separated domain walls in an otherwise ordered state is appropriate. In d = 1 such "kinks" 
destroy long-range-order at any non-zero T. The same crossovers also apply to the quantum ro- 
tor model of Section ^ the "ordered" regime is now one where spin-wave fluctuations about an 
ordered ground state dominate, and these fluctuations destroy long-range-order at any non-zero T 
for d = 2. 
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FIG. 5. Scaling function for the imaginary part of the dynamic susceptibiHty in the quantum 
critical region C of the transverse field Ising model in d = 1. The susceptibility x = Z^^/T"^/^ and 
k = ck/T, lJ = Lo/T. 
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FIG. 6. Scaling function for the imaginary part of the local dynamic susceptibility in the 
quantum critical region C of the transverse field Ising model in d = 1. The susceptibility Xl 
related to by (p5[), uJ = to/T, and notice the logarithmic scale on the horizontal axis. 
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FIG. 7. Generic phase diagram of a d-dimensional system displaying a second-order quantum 
transition at g = Qc and T = 0. The boundary of the shaded region represents a crossover and not 
a phase transition. See the text for more information. 
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